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ABSTRACT 

Although redshift-space distortions only affect inferred distances and not angles, they still 
distort the projected angular clustering of galaxy samples selected using redshift dependent 
quantities. From an Eulerian view-point, this effect is caused by the apparent movement of 
galaxies into or out of the sample. From a Lagrangian view-point, we find that projecting the 
redshift-space overdensity field over a finite radial distance does not remove all the anisotropic 
distortions. We investigate this effect, showing that it strongly boosts the amplitude of clus- 
tering for narrow samples and can also reduce the significance of baryonic features in the 
correlation function. We argue that the effect can be mitigated by binning in apparent galaxy 
pair-centre rather than galaxy position, and applying an upper limit to the radial galaxy sepa- 
ration. We demonstrate this approach, contrasting against standard top-hat binning in galaxy 
distance, using sub-samples taken from the Hubble Volume simulations. Using a simple model 
for the radial distribution expected for galaxies from a survey such as the Dark Energy Sur- 
vey (DES), we show that this binning scheme will simplify analyses that will measure baryon 
acoustic oscillations within such galaxy samples. Comparing results from different binning 
schemes has the potential to provide measurements of the amplitude of the redshift-space dis- 
tortions. Our analysis is relevant for other photometric redshift surveys, including those made 
by the Panoramic Survey Telescope & Rapid Response System (Pan-Starrs) and the Large 
Synoptic Survey Telescope (LSST). 
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1 INTRODUCTION 

The late-time acceleration of the expansion of the Universe has 
been one of the m ost exciting cosmological discoveries in recent 
years jRiess et al.|[T998 : Perlmutte r et al., 1999 ). Understanding the 
nature of this acceleration is one of the main challenges facing cos- 
mologists. One of the key observational methods that will be used 
to help meet this challenge involves using Baryonic Acoustic Os- 
cillations (BAO) in the 2-point galaxy clustering signal as a stan- 
dard ruler to make precise measurements of cosmological expan- 
sion. The acoustic signature has now been convincingly detected 



jPercival et alj200ll : ICoie et alj2005 *:'Eisenstein et al.'2005h using 
the 2dF Galaxy Redshift Survey (2dFGRS; CoUess et al. 2003) and 
the Sloan Digital Sky Survey (SDSS: lYork et al.ll2000l) . The de- 
tection has subsequently been refined using more data and better 
techniques, and is now producing intere st ing constraints o n cos- 
mological models (Percival et al. 2007a''bl: lOaztanaga et al.ll200^ : 
[Sanchez et al. 2009; Percival et al,2009.) . 

Some of the next generation of sky surveys, including the 
Dark Energy Survey (DES www.darkenergysurvey.org), 
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the Panoramic Survey Telescope and Rapid Response System 
(PanStarrs pan-starrs.ifa.hawaii.edu), and the Large 
Synoptic Survey Telescope (LSST www.lsst.org), will use 
photometric techniques to estimate galaxy redshifts, rather than 
more precise estimates from spectroscopic emission lines. The 
larger uncertainties on galaxy redshifts induce errors on inferred 
distances in the radial direction. The amplitude of the power 
spectrum and correlation function is reduced in the radial direc- 
tion by this smoothing, removing information. In this scenario, 
where little information remains from fluctuations in the radial 
direction, it makes sense to use the projected 2-pt functions in 
photo metric-redshift slices as t he statistics to com pare with mod- 
els tPadmanabhan et alj|2007l : iBlake etal]|2007h . The projection 
does not completely remove problems caused by inferring distances 
from velocity data (i.e. working in redshift-space). 

The distribution of galaxies that we observe in sky surveys, 
where we measure radial distances from spectroscopic or photo- 
metric redshifts, is not a true 3D picture. We observe an appar- 
ent clustering pattern in redshift-space, which is systematically dif- 
ferent from the true distribution in real-space because redshifts of 
galaxies are altered from their Hubble flow values by peculiar ve- 
locities. For example, on large scales, the infall of galaxies onto 
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Figure 1. The radial distribution of galaxies selected in a bin of width 
400 Mpc, calculated using photometric redshifts to estimate distances 
(solid line). This is compared against the distribution of true distances 
to these galaxies (dashed line) assuming a photometric redshift error of 
CTz = 0.03(1 + z). If the photometric redshifts of different galaxies are 
independent, then the expected projected conelation function of the photo-z 
selected sample, and a sample selected applying the dashed line as a selec- 
tion function based on the true distances, are the same. 



collapsed objects leads to an apparent enhancement of clustering 
in the radia l direction as galaxies are pro jected along their veloc- 
ity vectors ( lKaiseilll98l lHamilton|[T998h . When we infer galaxy 
distances assuming that the total velocity relative to the observer 
comes from the Hubble expansion flow, the result is that we see a 
distorted (redshift-space) density field. 

For angular measurements, these redshift-space distortions 
can alter the angular clustering in a redshift slice because the dis- 
tortions are correlated across the direction of projection. Although 
redshift-space distortions are sub-dominant compared with photo- 
metric redshift uncertainties, they give rise to a systematic effect, 
which nee ds to be included when photometric redshift surveys are 
analysed jPadmanabhan et alj|2007l : iBlake 67^120071) . This can 
complicate the analysis as the size of the redshift-space distor- 
tions, and therefore of this effect, is dependent on the cosmological 
model. Consequently, for every model to be tested against the data, 
we need to make a revised estimate of the redshift-space effect. 

In this paper, we consider the simplified problem in the plane- 
parallel approximation, and only consider linear redshift-space dis- 
tortions. Both photometric redshift errors and the random motion of 
galaxies in clusters provide an additional convolution of the over- 
density field along the radial direction. While these effects need 
to be corrected in any analysis, the required correction is easily 
modelled and can be separated from the linear redshift-space dis- 
tortion effects. For a measurement of the projected clustering, in- 
cluding such effects is equivalent to simply broadening the radial 
window function with which the galaxies were selected. This is 
demonstrated in Fig. [T] A top-hat bin in photometric redshift gives 
the same expected projected correlation function as simply apply- 
ing the convolved version of the bin as a selection function for the 
true distances. As we have to include a window function anyway, 
we simply assume in this paper that this window already includes 



the effects of both photometric redshift errors and the random mo- 
tion of galaxies in clusters. In the following analysis, we therefore 
assume that there are no redshift errors without loss of generality. 

The layout of our paper is as follows. In Section[2]we analyse 
the projected overdensity field and redshift-space effects upon it, 
both analytically (Section [2.1| & |23t and using Monte-Carlo sim- 
ulations ("Section [2.4t . We then consider how the recovered corre- 
lation function depends on galaxy binning (Section[3). Mock cata- 
logues drawn from the Hubble Volume simulation are constructed 
and analysed in Section |4] in order to validate this analytic work. 
We incorporate hybrid selection functions based on both real and 
redshift-space boundaries into the analysis in Section |5] In Sec- 
tion |6] we consider a non-uniform redshift distribution similar to 
that of future sky survey DES, and the realistic implementation of 
our work is discussed in Section|7] 



2 PROJECTED 2-POINT STATISTICS OF THE 
OVERDENSITY FIELD 

2.1 Correlation Function 

In order to simplify the problem, we assume that the clustering 
strength does not change across the samples under consideration 
and make the plane-parallel (distant observer) approximation, with 
redshift-space distortions along the z-axis of a Cartesian basis. In 
the absence of redshift distortions the projected correlation function 
is given by: 



(<5p(rp)5p(rp)), 



(1) 
(2) 



where d{rz,r'^,dp) = a/ {r^ — r^)^ + dp, and subscripts x, y 
and 2 denote the direction along each Cartesian axis, and p de- 
notes projected quantities p = xy. (j>{rz) is the radial galaxy se- 
lection function, normalised such that J dr^ fpirz) = 1, and 
^(d) = (5(r)5(r')), where S{r) is the overdensity of galaxies at 
real-space position r. Throughout our paper we use r to describe a 
galaxy position and d to describe the distance between two galax- 
ies, so, for example, is the position of a galaxy along the 2-axis, 
while dp is amplitude of the separation between two galaxies when 
projected into the x, y-plane. 

In reality, our radial position is determined via a redshift. In 
this case, Eq. ^ must be altered to 



Cp(dp) = ((5p(sp)(5p(sp)). 



(3) 



The weighted, projected overdensity field Sp{rp) can now be writ- 
ten 



l + <5p(rp) 



dsz (t>{sz)[l + S{s)], 



(4) 



where s = (rp, Sz) is the redshift-space position of each galaxy 
and <l>isz) gives the gala xy selection fu nction along the line of sight 
corresponding to Sz (e.g. fPeebledl 1980 1. 

The difference between the projection in redshift-space and 
real-space is shown schematically in Fig. [2] An edge to a window 
function (or a contour of constant galaxy density) that is straight in 
redshift-space is systematically distorted in real-space. The edge of 
the bin is itself clustered with a non-negligible projected correlation 
function, i.e. the real-space boundary has a correlation function that 
depends on rp . The inclusion or exclusion of galaxies is balanced 
in terms of the 3D correlation function within the boundary; while 
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Figure 2. Schematic showing the boundary of a region selected in redshift- 
space (solid line) compared with the boundary of the same region in real- 
space (dashed line). The boundary is distorted in real-space around an over- 
density and an underdensity. The positions of two galaxies whose apparent 
motion crosses the boundary are shown in redshift-space (solid circles) and 
in real-space (dashed circles). Note that, in this simplified picture where the 
under and overdensities have the same amplitude, the galaxy pair lost and 
the galaxy pair gained would contribute the same amount to the 3D real- 
space correlation function, following the dashed boundary. However, the 
projected clustering is different because we do not know the shape of the 
dashed line, and instead assume that the projection length is the same for 
all Tx^Ty. It is the 2D clustering strength of the boundary, and its correla- 
tion with the density field that is important, rather than the loss or gain of 
particular galaxy pairs. 



we lose voids, we gain clusters and these give the same clustering 
signal. However, we assume that the projected field has a constant 
projection length, and this implies that the underdensity of the void 
will become larger (since we include less of the galaxies) and the 
overdensity of the cluster becomes larger (since we will include 
more of its galaxies). Thus the overall clustering signal becomes 
stronger. 

The apparent shift in galaxy positions caused by moving from 
real to redshift space ( sz — rz) can be tr eated by Taylor expanding 
the selection function jFisheretal.ll993h . which gives to first order 



{Vz) + ■ 



drz 



-{Sz " Tz). 



(5) 



We consider this to be an Eulerian picture as it is based on apparent 
galaxy motions. We can write 



drz 



{rz)5(r) + [sz - Tz) 



d<t>irz) 



drz 



(6) 



to first order in S{r). Following linear theory, (sz — r^) can be 
written as a function of the overdensity field. 



{Sz 



(7) 



where /? = f/b, with / being the logarithmic derivative of the 
linear growth rate with respect to the logarithm of the scale factor, 
and b the galaxy bias. We therefore have that 



drz 



drz drz 



S{r) 



(8) 



If we think of (j)(sz) as setting up boundaries in Sz, then substitut- 
ing Eq. ^ into Eq. (O shows that we can expect coherent appar- 
ent galaxy motion across these boundaries. Correlations between 
galaxies moved into the sample by the redshift-space distortions, 
and those already within the sample, give rise to cross terms from 
the two terms in Eq. The second term in Eq. {8} also adds a 
component to the projected correlation function from the coherence 
of the velocities at different points on the boundary. We see that, 
even with constant (t>{sz) within a fixed interval, redshift-space dis- 
tortions can still affect the correlation function of the volume within 
the sample due to the motion of galaxies across the boundary. Mod- 
elling the effect of redshift-space distortions based on predicting 
galaxy motions (e.g. Regos & Szalay 1995) is difficult because we 
need to correlate multiple points on the boundary and internal lo- 
cations within the bin. 

In addition to the Eulerian picture given by Eq. {Sj, we can 
also consider a Lagrangian picture based on the redshift-space over- 
density field that we wish to project. Following this equivalent pic- 
ture, we can work directly with redshift-space overdensities using 
Eq. (H, 



Cpidp) = dszds'z(j){sz)4>{s'z)(.'' [d{sz,s'z,dp)\ 



(9) 



In the plane-parallel approximation, we can use the r edshift-space 
correlation function of equation 5 of lHamiltonl ( Il992h as input into 
the projection equation. 



r(d) = eo(d)Po(M) +6(d)P2(M) +C4(d)P4(M), 
where 

^o(d) = {b' + lbf + \f)ad), 



3 ■' 5- 



6(d) = 

C4(d) = ^.fm + ^^'{d)-;-^"{d)], 



Pi are the standard Legendre polynomials, and 



C'(d) 



r(d) 



3d" 



5d" 



i{d'){d'fdd', 



i{d'){d'fdd'. 



(10) 

(11) 

(12) 
(13) 

(14) 
(15) 



6 is the large-scale bias of the galaxy population being consid- 
ered, / is the standard dimensionless linear growth rate, ^ is the 
3-dimensional real-space correlation function, and ^ is the cosine 
of the angle between the separation along the line of sight and the 
transverse separation, ^. = \sz — s'z\/ d. One strong advantage of 
the Lagrangian framework is that it is straightforward to determine 
the projected correlation function, even when the galaxy selection 
function is discontinuous. This allows simple comparison between 
the results one expects to obtain with and without redshift-space 
distortions. 



2.2 The Limber approximation 

For pairs of galaxies, we can define the mean rriz = (rz + r'z)/2 
and separation along the z-axis dz = rz ~ r'z- For a survey whose 
depth is larger than the correlation length, and with a slowly vary- 
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ing selection function, so that ~ ^(r^) ~ (f){mz), Eq. ([2J 

reduces to tlie Limber equation in real-space (Sz ~ = 0) 



dd; 



+ di 



(16) 



We see that, for the Limber approximation, <;!) is a function of 
alone, and the integrals over drriz and ddz in Eq. U7b are separable. 
In redshift-space, a similar reduction of Eq. l|9} gives 



dm. 



dd. 



(17) 



if we expand redshift-space distortions in [sz — r^), or 

/ + 00 /' + 00 

dm, <^"(m,) / dd, r + - (18) 

oc ^ — oo 

in the Lagrangian picture. Because no galaxies are lost or gained 
moving from real-space to redshift-space, the result of the integral 
over dz is the same in real or redshift space, so we see that in this 
approximation there are no redshift-space effects. But, as we show 
later, this picture is too simplistic to be applied to the analysis of 
future data sets. 

2.3 Power Spectra 

InlPadmanabh an et al.l ( |2007|) . the projection of the 2-pt clustering 
was analysed through the power spectrum. We now consider such 
an approach in the plane-parallel approximation and for a Cartesian 
basis. Taking the Fourier transform of 5{s) in Eq. l|4) gives 



5,(rp) = j dsz 4>{sz) J ALj(k)e-"'-. 
We now define a window function 

\ —ikzSz 



W{kz) = / ds 



(19) 



(20) 



and use statistical isotropy and homogeneity within the definition 
of the power spectrum (5(k)5*(k')) = P(fc)5D(k — k'), where 
5d is the Dirac delta function. We assume that the power spectrum 
does not evolve over the volume covered by the windo\4j. Taking 
the 2-point function of the projected overdensity (Eq.|19b gives 



Cp(dp) = {5p(rp)(S(rp)) 



-W\kz)P{k)e 



-ikp.(rp-r' ) 



(21) 

(22) 



(27r)3 

The projected overdensity can be written in terms of a 2D power 
spectrum Pp{kp), 



^p(dp) 



(23) 



(24) 



If we compare Eqns. J22b & J23t , we see that 

Note that the power P(k) depends on the amplitude of the full 3- 
dimensional wavevector, and so is dependent on fcp. 

Using Eq. (|5} to include redshift-space distortions, the window 
W{kz) has an extra term, 

^ This is true if analysing a single time slice from a simulation 



W(kz 



dr. 



dr. 



g->fc,r,^ (25) 



In Fourier space, (s, — r,) — ~/3(fc^/fc'^)(5(r), so we can ex- 
pand (5(s) to 1st order in 5(r), leaving a new window function for 
Eq. l|24j 

■ ' d4,(r 



W{kz 



drz 



dr. 



e " (26) 



If we drop the plane-parallel approximation and expand in 
Spherical Harmonics, the standard result ( Peebles. 1973i) is 



(|a!m| > 



where 



27r2 



dk k^P[k)W'^(k), 



W{k)= I dr ^{r)Mkr) + ^^jiikr) 



(27) 



(28) 



Here the I dependence is contained within W{k), while in Eq. l l24b . 
it was the power that depended on kp. Eq. i24\ could have been 
rewritten by changing the variable of the convolution integral to k 
to match. 



2.4 Monte-Carlo simulations of the projection effect 

In order to test the projection formulae presented in Sections 12.11 
& 12.31 without redshift-space distortions, we have used Monte- 
Carlo realisations of ^- function real-space co rrelation functions in 
a similar vein to that o f ISimpson et"ai] j2009l) . We work in a plane 
parallel approximation throughout and construct a real-space 3D 
(5-function correlation function at an arbitrary location do such that 



C(d) = fo(d-do)Co, 



(29) 



where Sd is the standard Dirac delta function. We do this by intro- 
ducing a pre-determined excess of data pairs at the location do. 
The number of excess pairs we introduce depends on the value 
of ^(do) we require and is determined using the natural estima- 
tor ^ = D/i? — 1. For example, if we have a uniform distribution 
of data and random pairs with 100, 000 pairs per bin of separation, 
we would require an excess of 10, 000 data pairs at the location 
do for ^(do) = 0.1. In doing this we create an unnormalised 3D 
5-function correlation function. 

Changing the variables in the inner integral of Eq. il6\ to be a 
function of 3D pair separation d gives 

Udp)= f f dm.dd (l>\mz) (30) 
J Jv ^Jd^ ~ dl 

and is simplified for the 5-function case such that 

^p(dp) = ^ i dniz 0^(m,)go f° . (31) 
J ^dl - dl 

The factor l/vrdo accounts for the fact that the (5-function real- 
space correlation function was unnormalised. By introducing a ra- 
dial window, we are preferentially selecting pairs of galaxies from 
the sample. A further volume reduction normalisation is required in 
Eq. to account for this. The excess probability of finding two 
galaxies in areas SAi and 5A2 with a 2D projected separation dp is 
the sum of all the probabilities of finding two galaxies in volumes 
SVi and 5Vj along the radial axis at all 3D separations d. That is. 



1+Cp(dp) 



nv 



2^2^[i + ad.,)]SV,SV, .(32) 
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Figure 3. Projected con'elation functions calculated for a 3-dimensional 
(5-function correlation function centred on d-j,D = 90 /i^^ Mpc, with no 
radial window (solid symbols) and with a top-hat window in radial distri- 
bution of width 100 Mpc (open symbols). Models calculated using 
Eq. j3U are shown by the solid Unes. 
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Figure 4. A schematic representation showing how galaxy pairs are se- 
lected using top-hat and pair-centre binning schemes. Using a top-hat bin- 
ning scheme, where galaxy pairs are selected according to the position of 
each individual galaxy, pair AB would be placed in redshift bin 2, whereas 
pair CD would not be placed in any bin, and would simply not be counted 
in an analysis. In contrast, the pair-centre binning scheme would place both 
pairs in bin 2. 



In Fig.|3]we show the clustering expected for a projection of 
a density field created from a 5i5 -function 3D correlation function 
in the case where there is no window function (solid symbols) and 
for a window function of width 100 h^^ Mpc (open symbols). The 
excess of pairs that exists at a single scale in 3D is projected onto 
a range of scales, up to and including this scale, in 2D. The projec- 
tion window leads to a damping of power on all scales. This effect 
depends upon the window size; as we move to smaller projection 
windows the effect of the projection is decreased and ^215 — > £,3d- 
The projection of a more general density field, where there is clus- 
tering on a range of scales, can be considered as the linear combi- 
nation of the projections of a series of (Ju-function 3D correlation 
functions. The trends observed in this analysis will help us to inter- 
pret the behaviour of the projected correlation function in the more 
general situation analysed in later sections. 



3 BINNING GALAXY SAMPLES 

Future surveys will automatically have a standard selection func- 
tion caused by the changing cosmological volume, the number 
density of galaxies as a function of redshift, and selection effects 
such as a magnitude limit below which we cannot observe galax- 
ies or obtain accurate photometric redshifts. In addition to this dis- 
tribution we will wish to bin galaxies based on their photomet- 
ric redshifts in order to analyse the evolution of galaxy proper- 
ties and/or cosmology across the sample. We now consider how 
the way in which this sub-division is applied affects the importance 
of redshift-space distortions. 

One simple approach would be to bin galaxy positions in red- 
shift, equivalent to a top-hat binning. Such galaxy selection means 
that galaxy pairs, where galaxies lie in different bins, are not in- 
cluded in the estimate of the correlation function. This exclusion of 
pairs leads to the observed difference between the projected real- 
space and redshift-space correlation function, as described in Sec- 



tion l2.1l An alternative to this approach, considered here, would be 
to bin galaxy pairs rather than individual galaxies. 

A simple argument shows that in an ideal situation, applying 
a binning based on the centre of galaxy pairs in the radial direc- 
tion, which hereafter we refer to as pair-centre binning, can com- 
pletely remove the effect of redshift-space distortions while retain- 
ing information about the evolution of the correlation function. A 
schematic representation of this binning scheme is shown in Fig.|4] 
Consider the galaxy pair defined by galaxies A and B: the posi- 
tions of both galaxies and their pair-centre are within redshift slice 
2. This pair would therefore be included in analyses conducted on 
this slice in both top-hat and pair-centre binning schemes. The po- 
sitions of galaxies C and D span two separate redshift slices and 
therefore the pair they define would not be included in an analysis 
of either slice 2 or 3 when using a top-hat binning scheme. How- 
ever, this pair would be included in an analysis of slice 2 when 
using the pair-centre binning scheme. This schematic demonstrates 
both the pair-centre binning scheme and the fact that such a scheme 
includes all pairs within an analysis. 

Suppose that we have a clustered distribution of D galaxy 
pairs of separation r with a uniform sampling function along the 
z-axis in a large volume that would contain R pairs if galaxies were 
randomly distributed. Because of the large volume assumption, we 
can assume that boundary effects for this sample are negligible. 
Therefore, redshift-space distortions have no effect for the full cat- 
alogue for which our estimate of ^p{dp) is (,p{d,p) = D/R — 1. 
Now suppose the sample is split into n sub-samples, based on the 
redshift-space positions of the centres of the pairs within equal vol- 
umes, chosen independently of the observed galaxy distribution. 
Then all pairs are still counted in some bin; none are lost or gained 
as opposed to galaxy based selection functions. For the subsam- 
ples, {D') = D/n, R' = R/n, and {£,p{d,p)) is unchanged from 
the value for the full sample. This is true regardless of bin size. The 
key difference here, compared with considering a set of bins based 
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Figure 5. Top panel: the normalised radial distribution of galaxies (solid 
line) and pair-centres (dashed Hne) for the distribution of galaxies in a top- 
hat bin of width 100 h^^ Mpc. These are compared with the distributions 
of galaxies (dot-dash line) and pair-centres (dotted line) for galaxies whose 
pair-centre is within a 50 Mpc bin, and with < 100 Mpc. 
Bottom panel: comparison of the radial pair separations (ciz), between top- 
hat (solid line) and pair-centre (dashed line) binning. 



on galaxy selection, is that no pairs are left out, so the expected 
correlation function has to be the same for all bins. 

For a sample where we do not know the true distance to each 
galaxy, but instead rely on photometric redshifts, binning based 
on apparent pair centre will also remove redshift-space distortions. 
The above argument based on pair conservation will also hold in 
this situation. 

The radial distributions of galaxies and pair-centres along with 
the distribution of radial pair separations for top-hat and pair-centre 
binning schemes are compared in Fig. [5] For measuring the radial 
evolution of clustering through binned projected correlation func- 
tion measurements, there is no obvious advantage to either scheme. 
This is particularly true when photometric redshifts are used to es- 
timate radial positions, as it is then impossible to select galaxies 
from non-overlapping radial bins (see Section|7]l. 

We therefore see that we can add boundaries based on pair- 
centres and analyse projected clustering in bins without being af- 
fected by redshift-space distortions. However, there are two prob- 
lems with applying this approach in practise: 

(i) galaxy pairs of wide separation now have to be included, 



(ii) galaxy surveys typically have flux limited boundaries, which 
cause redshift dependent effects that cannot be removed by any 
binning. However, this effect can be removed by fc-correcting the 
observed luminosities and cutting the sample at a more stringent 
fc-corrected luminosity limit. We now investigate this further. 

3.1 Flux-limited Selection Functions 

Peculiar velocities can directly influence galaxy brightness through 
relativistic beaming, but such effects are small for typical galaxy 
peculiar velocities. Redshift distortions would additionally change 
the apparent magnitudes through the fc-correction, potentially caus- 
ing galaxies to either enter or exit flux-limited samples. The change 
in apparent magnitude will correlate with bulk-flow motions and 
thus the boundary of the survey in real-space will fluctuate in a 
manner analogous to that described in Fig. |2] In this situation, the 
amplitude of the effect and whether it enhances or reduces the real- 
space clustering signal will depend on galaxy type and the band 
used for detection, but for a homogeneous sample of galaxies (e.g. 
Luminous Red Galaxies) one would expect that this effect will be 
significant. 

This redshift-space effect is simple to remove - fc-corrections 
derived by fitting to galaxy spectra will correct for spectral shifts 
caused by both the Hubble flow and any peculiar velocities. It 
therefore makes sense to select galaxy samples after applying the 
fc-correction, and cutting back from survey boundaries based on 
apparent magnitude, until no galaxies outside the original sample 
would be expected to pass the revised boundary. This is not as oner- 
ous as it sounds as one has to do this to create true volume-limited 
catalogues. 

Given purely photometric data, fc-corrections can only be es- 
timated given a photometric redshift and spectral-type fit, and are 
therefore unreliable for individual galaxies. For this reason, and the 
fact that cutting back from the survey boundary removes a large 
amount of data, fc-correction s have not always been applied to ap- 
parent magnitude limits (e.g. lRoss & B runner 2009 select galaxies 
with de-reddened r < 21 for their parent sample). We therefore 
now consider the amplitude of the effect. One can express the fluc- 
tuation in magnitude, 5m, as 

5m = Akcorr/dz5z (33) 

where 5z is the magnitude of the redshift distortion. This will cause 
fluctuations in the effective depth of the survey such that 

DM{zeff) ~ DM{z) = 5™. (34) 

where DM{z) is the distance modulus, Zejf is the effective depth 
and 2 would be the predicted depth. The SDSS DRV photomet- 
ric redshift table includes estimated r-band fc-corrections for ev- 
ery galaxy. Studying galaxies with type- value equal to (the most 
early-type), one can determine that dkcorr/dz ^ 3.3 at z — 0.4. 
For an arbitrary 5z, this dkcorr/dz yields z^ff ~ z = 0.5Sz. 
For example, assuming bulk flows have a velocity ~ lO^km/s — 
thereby imparting redshift distortions at the ~ 1% level (5z — 
0.004) — they impart coherent fluctuations in apparent magnitude 
equivalent to 0.013 magnitudes (in the r-band). At z = 0.4, these 
fluctuations in magnitude imply a change in the survey depth of 
Zeff — z = 0.002 (0.5%). Thus, the redshift distortions caused by 
selecting a flux-limited sample of galaxies can be as large as 50% 
of those caused by selecting a sample in redshift. Therefore, even 
for a flux limited selection function, redshift distortions may be im- 
portant. The size of the effect depends on the slope of kcorr{z), and 
one can minimise the effect by carefully choosing the band used for 
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selection and the type(s) of galaxies included in the sample. (One 
can envision cases where slope of the average fc-correction is zero, 
thus removing any effect even before applying fc-corrections.) 



4 ANALYSIS OF HUBBLE VOLUME SIMULATIONS 

In order to test the effect of redshift-space distortions on the pro- 
jected correlation function for a realistic non-linear distribution of 
galaxies, we have analysed results fro m the ACDM Hubble Volume 
(HV) simulations l lEvrard et alj2002h . The ACDM HV simulation, 
covermg a (3000 /i-^Mpc)^ box, assumes a cosmological model 
with n„ = 0.3, Q.CDM = 0.25, fib = 0.05, = 0.7, h = 70, 
as = 0.9, &n, = l. 

We make a number of simplifications in order to help with 
the calculation of projected real-space and redshift-space correla- 
tion functions. For each sample to be analysed, along the two non- 
projection axes, we use the periodic nature of the numerical simu- 
lation to eliminate boundaries. This means that we can confidently 
use the natural estimator ^ + 1 = D/ R, where the expected number 
of galaxy pairs in the absence of clustering R can be calculated ana- 
lytically. We also do not introduce a galaxy-bias model, and assume 
that galaxies Poisson sample the matter particles. The inclusion of 
such a model would not alter the conclusions of this work. 

We start by applying a top-hat selection function to the galaxy 
positions, calculating projected correlation functions for window 
widths 5Qh^^ Mpc and 100 /i^^ Mpc in real and redshift space. 
Fig. [6] shows the correlation function after reducing noise by 
averaging over 30 samples (100 Mpc bins) or 60 samples 
(50 Mpc) bins. In real-space the projected correlation function 
tends towards the 3D correlation function at large scales, as ex- 
pected. In line with the analysis presented in Section [2!4l the scale 
at which becomes ~ ^3_d is larger for the 100 h^^ Mpc bin. For 
each bin size, the inclusion of redshift-space distortions clearly has 
a strong effect and this effect grows dramatically as the scale gets 
larger. Notably, it is larger even than the effect of redshift-space 
distortions on the 3D spherically averaged correlation function (or 
power spectrum). The effect is enhanced in the narrower projection 
window. As well as increasing the amplitude of the projected cor- 
relation function, we see that redshift-space distortions also act to 
wash out the baryon acoustic oscillation signal. 

Selecting galaxy pairs solely based on the position of their 
pair-centre removes the effect of redshift-space distortions. To see 
this, suppose we split along the projection axis into A'^ slices, and 
average the DD counts over all slices. Then the average is inde- 
pendent of as all pairs are counted however many bins are se- 
lected. In addition, the periodic nature of the simulation means that 
no pairs are gained or lost between real-space and redshift-space: 
we always count all pairs of galaxies, so there will be no change 
in the measured correlation function. As explained in Section [3T1 
if we select based on an apparent magnitude limit, we can remove 
redshift distortions by applying a more stringent magnitude limit 
based on fc-corrected luminosities. Here we have to cut the lumi- 
nosity limit back to make sure that the new sample is complete, 
in that it contains all of the possible galaxies. However, there is a 
further practical problem in that including galaxy pairs with wide 
radial separation might complicate the modelling of cosmological 
evolution required to fit the correlation function. Consequently, it 
might be difficult to analyse the measured correlations function for 
a pair-centre binned sample in practice. 

We therefore introduce a constrained pair-centre binning 
scheme that includes an upper limit on the pair separation along the 
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Figure 7. The expected ratio of the projected correlation functions in 
redshift-space and in real-space, averaged for "angulai'" sepai'ations be- 
tween 40 h~ ^ Mpc and 80 /i~ ^ Mpc, as a function of bin width. The solid 
line show the difference as a function of the width of the top-hat window. 
The dashed line show the result for constrained pair-centre binning as a 
function of an additional constraint placed on the radial galaxy separation. 
We have plotted results (and therefore matched filters) as a function of the 
mean radial galaxy separation. 



projection axis, in addition to pair-centre binning. This is equiva- 
lent to locating each galaxy included in the analysis in the centre 
of a top-hat bin. We should expect that the effect of redshift-space 
distortions will be reduced compared with binning galaxy distri- 
butions in a top-hat with the same width, as boundaries will only 
affect galaxy pairs with the maximum radial separation, whereas 
for top-hat bins they affect galaxy pairs with a range of radial sep- 
arations (see Fig. |5j. Results calculated using this binning scheme 
are shown in Fig. |6] Here we see that the effect of redshift-space 
distortions is reduced, especially for the larger \dz \ limit. 

In order to investigate the effect of different binning schemes 
further, Fig.|7]shows a comparison on the large-scale redshift-space 
and real-space correlation function amplitude. These are averaged 
for galaxy separations between 40 Mpc and 80 Mpc. We 
have plotted these as a function of average radial galaxy separa- 
tion, in order to compare filters in an unbiased way. We clearly see 
that, when binning radially using the constrained pair-centre bin- 
ning scheme, the effect of redshift-space distortions is significantly 
reduced. 

The relative importance of redshift-space distortions depends 
on the average galaxy bias of the populations being considered; 
there is a balance between the impacts of h and / in Eq. l llOt . In 
order to demonstrate this. Fig. [8] shows that the relative effect of 
redshift-space distortions decreases as the bias of the galaxy sample 
analysed increases. This explains why the effect of redshift-space 
distortions was reduced in the work of Baldauf et al. (2009). 

In this section, we have considered the cases of a top-hat or 
pair-centre galaxy selection. We have argued that while, in prin- 
ciple, pair-centre binning removes the effects of redshift-space dis- 
tortions provided fc-corrections are included when magnitude limits 
are applied, there are good reasons to remove galaxies of wide sepa- 
ration if we are to measure the evolution in the correlation function. 
Therefore, in Section |6] we will test how these binning schemes 
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Figure 6. Top row: Correlation functions calculated from Hubble Volume data in radial bins of width 50h~^ Mpc and iQOh~^ Mpc. Solid symbols are 
plotted where the correlation function is positive, while open symbols show where the correlation function is negative. Bottom row: Correlation functions 
calculated from Hubble Volume data for galaxy pairs selected based on constrained pair-centre binning scheme with both their pair-centres and radial separation 
less than 50 /i^^ Mpc or 100 /i^^ Mpc. Icr error bars are plotted in both cases, assuming that the slices analysed draw coiTelation functions from a Gaussian 
distribution. The dotted line gives the 3D HV correlation function, liuv^ (plotted assuming r = Vp) as measured from the simulation. The solid lines denote 
the projected correlation function one expects in real-space with ^hv as the 3D correlation function (lower curves in each panel), and in redshift-space using 
Eqns. (9) & (To) to estimate the 3D redshift-space coiTelation function from iiuv (upper curves in each panel). 



work when the background galaxy distribution has a more realis- 
tic radial distribution, similar to that expected for a survey like the 
Dark Energy Survey. Before we can do that, we need to consider 
the case where we have a boundary that consists of a mix of real- 
and redshift-space constraints. 



5 DEALING WITH HYBRID SELECTION FUNCTIONS 

In practice, the radial selection function will be dependent on both 
observational constraints such as the limiting apparent magnitude 
of the survey, and additional binning. One expects that the bound- 
ary based on observational constraints can be treated as a real-space 
boundary (though, this even, is not so simple; see Section 13. 11 1. 
Thus, when one applies a top-hat selection in redshift to an ob- 
served sample of galaxies, the resulting boundaries of the selec- 
tion function will include both real-space and redshift-space com- 
ponents. 

Fig.|9]shows a schematic representation of a top-hat selection 



in redshift made at positions s^i and a^2 along a non-uniform real- 
space radial selection function. It shows that we can split galaxies 
within this bin into three sub-samples, with different boundaries: 

• As (redshift-redshift): Selected with both boundaries in 
redshift-space. 

• Bh (redshift-real): Selected with one boundary in real-space 
and one boundary in redshift-space (/lyind-space). 

• Cr'- Selected with both boundaries in real-space. 

The real-space and redshift-space boundaries of Fig. [9] are repre- 
sented by solid and dotted lines respectively. Any auto-correlation 
of galaxies with this selection function will essentially be a 
weighted sum (based on the amplitude of the selection function) 
of the auto-correlations of galaxies within the individual subsam- 
ples and the cross-correlations of galaxies in different subsamples. 

In order to investigate the projected clustering of these dif- 
ferent subsamples, we have drawn samples of particles from the 
HV simulation (see Section |4l(, created in top-hat bins of width 
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Figure 8. The expected ratio of the projected correlation functions in 
redshift-space and in real-space, averaged for "angular" separations be- 
tween 40 /i^^ Mpc and 80 Mpc, as a function of galaxy bias, as- 
suming a ACDM cosmology with Qrn = 0.25. The solid line show the 
difference as a function of the width of the top-hat window. The dashed line 
show the result for constrained pair-centre binning as a function of the ad- 
ditional constraint placed on the radial galaxy separation. As in Fig. [7] we 
have plotted results (and therefore matched filters) as a function of the mean 
radial galaxy separation. 











As 





Figure 9. Schematic representation of an evolving real-space radial selec- 
tion function with populations As , i?h and C'r defined according to where 
a top-hat bin with redshift-space boundaries at s^i and 8^2 intersect the 
radial selection. Populations have boundaries in: As redshift-space, Bf^ 
hybrid-space and Cr real-space. 



Sample As has top-hat selection boundaries in 
redshift-space, sample has one real-space and one redshift- 
space boundary, while sample Cr has both boundaries in real- 
space. These samples cover the same region of the simulation. 

Fig. [To] shows the projected auto-correlation functions for 
these subsamples. The measured for the Cr and As samples are 
essentially the same as those shown in the top-right panel of Fig. (6] 
and just as before they return the expected real and redshift-space 
correlation functions calculated via Eqns. ^ & l IlOl l. However, the 
hybrid-space correlation function, of sub-sample Bh has an am- 
plitude that lies in-between those of the pure real and redshift-space 
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Figure 10. The average recovered auto-correlation function (solid circles) 
for galaxies from 90 samples drawn from the Hubble Volume simulation us- 
ing three different radial selections, each with top-hat width 100 h~^Mpc. 
These are compared against model correlation functions calculated for dif- 
ferent galaxy samples Eq. (36). The three radial selections are: 1) two 
real-space boundaries (lowest points), which best matches the model calcu- 
lated using the real-space correlation function, 2) two redshift space bound- 
aries (highest points), which best matches the model calculated using the 
redshift-space correlation-function and, 3) a real-space boundary on one 
side and a redshift space boundary on the other side (points in the middle), 
which best matches the model calculated using the geomefilc mean of the 
real- and redshift-space congelation functions. 



correlation functions. We find that we can effectively model 
by assuming the underlying 3D overdensity field has a correlation 
function given by 



(35) 



Note that we are using 5' to represent the real-space 3-dimensional 
correlation function. As can be seen in Fig.[To] this model is well- 
matched to the measured ^p. The justification for this model is that 
the multiplicative boost to the projected density fluctuations (R if 
we consider that ^ = D/R — 1) can be decomposed into multi- 
plicative contributions from each boundary. Following this model, 
we should find that the relative effect of redshift-space distortions 
on each population, and their cross-correlations, are simply pro- 
portional to the number of redshift-space boundaries present. If we 
choose galaxies from a sample with n £ {0, 1, 2} redshift-space 
boundaries, and another from a sample (possibly the same one) 
with m G {0, 1, 2} redshift-space boundaries, then expected corre- 
lation function is given by 



(36) 



where I = m + n. 

Fig. [TT] displays the cross-correlations between the our three 
HV subsamples. As expected, the model calculated using the ap- 
propriate 5'' from Eq. i36\ is the closest match to the measured 
cross-correlation in every case. All of the models do over-predict 
all three measurements at large scales, but we believe this is re- 
flective of the error associated with our measurements (one would 
expect it to be covariant between each sample as they all sample the 
same density field). It is possible that we are seeing effects caused 
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Figure 11. The average measured cross-correlation functions from 90 radial 
slices of width 100 h^^Mpc in real-space, redshift-space or a hybrid with 
one real-space and one redshift-space boundary, each containing 10'' galax- 
ies (solid circles). These are compared against the model of Eq. )36t 
(solid hnes), for different total numbers of redshift boundaries. The ampli- 
tude of both model and data correlation functions increase with increasing 
dependence on the redshift-space correlation function. 

by the coherence of the boundaries with each other that would be 
removed for wider bins, such as those we consider in Section|6] 

Given a hybrid selection function such as that shown in Fig.[9l 
we must split the sample into populations where we can assume 
simple boundary conditions for each. In fact, we can consider 
solving the projection equation (e.g. Eqns. [2] & |9] in real-space 
and redshift-space) by Monte-Carlo integration over pairs of radial 
galaxy locations. For each pair of locations we can determine the 
relative contributions from galaxies in each of the subsamples, and 
therefore construct a full model for the correlation function. 



6 IMPLICATIONS FOR FUTURE PHOTOMETRIC 
REDSHIFT SURVEYS 

A number of extremely wide angle imaging surveys are planned 
over the next few years: the Dark Energy Survey (DES), the 
Panoramic Survey Telescope & Rapid Response System (Pan- 
Starrs) and the Large Synoptic Survey Telescope (LSST). One goal 
of these surveys is to constrain the current acceleration of the Uni- 
verse. In general, one can hope to use such surveys to make four 
measurements of dark energy using complimentary techniques: 
cluster counting, BAG, weak lensing and supemovae. In this pa- 
per we consider BAG measurements. For these experiments, radial 
distances to galaxies will be estimated from photometric redshifts, 
so there will be little information in the radial direction on the scale 
of BAG. Consequently, analyses will tend to rely on making pro- 
jected galaxy clustering measurements in redshift slices that are 
sufficiently narrow to be able to reveal cosmological acceleration. 

In order to assess the effect of redshift-space distortions on 
such measurements, we now consider one of these surveys, DES, 
in more detail. The DES will use a 500 Mega-pixel camera on the 
Blanco 4-metre telescope in Chile to conduct a galaxy survey over 
a sky area of 5000 deg^ . Multi-band observations using g, r, i and z 




Figure 12. Approximate redshift distribution similar to that expected from 
the Dark Energy Survey. In order to use this distribution of galaxies to eas- 
ily measure cosmological acceleration using projected clustering measure- 
ments, this population will have to be subdivided or binned in redshift. 
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Figure 13. Top panel: normalised radial selection functions for top-hat 
slices of width 400 Mpc created from a DES-like distribution. Bot- 
tom panel: We also consider bins in radial galaxy pair-centre of the same 
width 400 h-^ Mpc. While we bin in distances derived from photomet- 
ric redshifts (solid lines), the true distiibution of radial galaxy distances is 
shown by the dashed lines. 
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filters will allow photometric redshifts to be obtained over a range 
0.2 < z < 1.4. The expected redshift distribution of the galaxies 
will be approximately 




after applying approximate survey depths to basic luminosity func- 
tions. This function is plotted in Fig. [12] This distribution of galax- 
ies will then be sub-divided into bins in order to assess the evo- 
lution of the BAO scale across the survey. As discussed above, 
measurements of the projected correlation function will be af- 
fected by redshift-space distortions, which will increase the signal 
strength and decrease the importance of BAO features. We now 
consider how the choice of binning methodology affects the impact 
of redshift-space distortions. 

We consider splitting this galaxy distribution into five redshift 
slices each of width 400/i~^Mpc for distances estimated from 
photometric redshifts, assumed to be Gaussian with = 0.03(1 + 
z). These bins cover radial distances of 500 — s- 2500 Mpc, re- 
lated to redshifts z — 0.15 to 2 = 1.06 (assuming a flat ACDM 
cosmology with Q,m ~ 0.25. The upper panel of Fig.[T3]shows the 
distributions of galaxies in these slices. The lower panel of Fig.ll3l 
shows the redshift distributions when we bin the galaxies based on 
the centre of the radial separation, calculated from the photometric 
redshifts. Because we are using photometric redshifts, there is no 
way to bin without leaving overlap in the true radial distributions. 
Consequently, the top-hat binning scheme does not provide an ob- 
vious advantage over other schemes in terms of analysing disjoint 
regions. 

In light of the discussion in Section lBTT] we consider both the 
case in which ijiu B5 is treated as a real-space boundary (results pre- 
sented in Fig. 114b. and the case in which it is treated as a redshift- 
space boundary (as may be the case when the slope of kcorr{z) is 
especially large; results presented in Fig.llSb. For the hybrid bound- 
ary, we employ the techniques described in Section[5]to determine 
the full form of the projection. When we treat <j!>_d_bs as a redshift- 
space boundary, we can simply use Eq. (|2l( & ([9} to determine 
in real and redshift-space. 

The left-hand panels of Figs. 1141 and 1151 show the expected 
projected correlation functions when a top-hat binning scheme is 
applied with width 400 Mpc. Even for this large bin width, 
in every radial bin there is a significant difference between the re- 
sult obtained using the redshift-space correlation function and the 
real-space correlation function. The difference is made clear by ob- 
serving the ratios between the two, displayed in the bottom panels. 
The ratios are slightly higher in the case where we treat (j>DES as 
a redshift-space boundary (Fig.llSt. and the difference between the 
two treatments is largest for the lowest redshift bin (which has its 
selection most affected by the overall DES selection). In every case, 
the ratio is significant 1.5) around 100 Mpc and the shape 
of the predicted and measurements differ substantially. 

The effects of redshift distortions are completely removed 
when a pair-centre binning scheme is employed and the cj>DES 
boundary is assumed to be real-space, as made clear in the middle 
panel of Fig. [141 Based on the discussion in Section[3] we can sim- 
ply use Eq. ([2} for both and thus their ratio is identically I. Even 
when the <j)DES boundary is assumed to be in redshift-space, as 
displayed in the middle panel of Fig. [15] the difference between the 
redshift-space and real-space model is considerably smaller than 



^ We thank the DES LSS working group for providing this approximation 



for the top-hat binning. The right-hand panels of Figs. 1 141 and 1 151 
show that even if one applies the constraint that the separation be- 
tween pairs be less than A00h~^ Mpc to be included in a con- 
strained pair-centre bin, redshift space distortions introduce a much 
smaller effect than for a top-hat binning scheme. 



7 DISCUSSION 

Redshift distortions produce a strong effect on projected clustering 
measurements — one that is far stronger than the redshift-space dis- 
tortion effect on the 3D clustering signal for galaxy samples with 
low bias and a narrow radial window. It is clear that redshift distor- 
tion effects must be included when modelling the projected galaxy 
clustering in redshift slices. 

If we consider the apparent motion of galaxies as we move 
from real- to redshift-space, then redshift-space distortions cause 
an apparent coherent motion of galaxies into and out of samples. 
This is true whether samples have sharp boundaries, or if the se- 
lection function changes more gradually with distance. In fact, we 
have argued that such motion does not in itself alter the projected 
correlation function — we would recover the real-space projected 
correlation function if we could correct for the movement of the 
boundary (i.e. allow for the depth of the survey to change with the 
distortions). However, this is not easy to do, although it is theo- 
retically possible and is an interesting alternative approach. The ef- 
fect of redshift-space distortions is due to the redshift-space bound- 
aries themselves having an angular clustering signal, and their cor- 
relation with the overdensity field. We can alternatively view the 
effect from a Lagrangian standpoint, where we have to consider 
that the projection does not remove redshift-space effects from the 
anisotropic correlation function. 

We have used Hubble Volume simulations to show that the 
projected correlation function can be modelled most easily by in- 
tegrating the redshift-space correlation function over the radial se- 
lection function. Galaxy selection will often be a mix of real and 
redshift-space constraints, and we have shown that this can be mod- 
elled by splitting the population into samples that can be considered 
to have top-hat windows in either real-space, redshift-space or a 
hybrid of the two. In the hybrid situation, the projected correlation 
function can be modelled using both the real-space and redshift- 
space correlation function over the radial selection function, and 
that more complicated selection functions can be effectively mod- 
elled in a similar manner. Prior to this publication, no-one has con- 
sidered how these hybrid selection functions affect the recovered 
projected clustering signal. 



7.1 Pair-Centre Binning 

We have presented a new measurement technique, pair-centre bin- 
ning, and shown that it minimises the effects of redshift space dis- 
tortions. In this new scheme, we only include galaxies where their 
apparent pair-centres lie within a given radial bin, whereas tradi- 
tional methods select pairs where both galaxies lie within the bin. 
The new scheme includes individual galaxies that lie outside the 
traditionally applied top-hat boundaries. This simple modification 
acts to reduce the effect of the coherent movement of galaxies be- 
tween slice boundaries on projected correlation function clustering 
analyses. It is important to note that this new technique does not 
prevent the movement of galaxies between slices; redshift-space 
distortions due to peculiar velocities will always exist in the radial 
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Figure 14. Top panels: Real-space (dashed lines) and redshift-space (solid lines) coiTelation functions predicted for the 5 radial bins drawn from the DES-like 
selection function, assuming it can be treated as a real-space boundary. Bottom panels: The ratio between the redshift-space and real-space projected correlation 
function. Here different line styles correspond to different bins: in the order of increasing redshift, they are solid, dashed, dot-dash, dotted, dot-dot-dash. From 
left to right: Top-hat bins of width 400 h^^ Mpc in the radial direction, pair-centre bins of width 400 h^^ Mpc, and constrained pair-centre bins of width 
400 h^^ Mpc, including an additional constraint on the radial separation of \dz \ < 400 Mpc. 
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Figure 15. As Fig. 1 141 only in this case we treat the DES selection function as a redshift-space boundary. 



direction. It simply makes sure that they do not produce a coherent 
effect on the measurements. 

There are two potential disadvantages of the pair-centre bin- 
ning scheme. One is the fact that the same galaxy may be included 
in multiple radial bins — thus introducing a correlation between 
radial bins. Another is the fact that such a scheme results in nec- 
essarily wider radial bins, which causes the clustering signal to be 
diluted. We do not feel that either is a large problem. Applying 
the more traditional top-hat binning scheme to photometric surveys 
necessarily results in overlapping radial bins (due to photometric 
redshift errors) and there will always be considerable covariance 
between radial bins selected with photometric redshifts — we do 



not think that pair-centre binning will make this problem consid- 
erably worse. The dilution effect can be mitigated by imposing a 
maximum separation between the pairs included in a pair-centre 
bin: we call this constrained pair-centre binning. As can be seen 
by comparing the middle and right-hand panels of Figs. [141 and[T5| 
imposing such a constraint increases the expected signal while not 
causing a significant change in the effects of redshift-space distor- 
tions. More detailed studies of these effects are warranted, but we 
are confident that the reduction in the redshift distortion effect we 
observe when utilising pair-centre binning will make this scheme 
considerably preferable to a top-hat binning scheme. 

Pair-centre binning completely removes the effect of redshift 



© 2008 RAS, MNRAS Q0Q,[TlfT3] 



Redshift-space distortions and projected clustering 1 3 



distortions when given a uniform galaxy distribution. Sucli perfect 
distributions do not exist — most galaxy samples selections are 
based on an apparent magnitude limit — and thus realistic radial 
distributions of galaxies are more complicated. However, we have 
argued that if galaxy samples selected based on an apparent magni- 
tude limit are cut back so that no galaxies fc-corrected galaxies are 
missing from the sample, then this does not matter: the boundaries 
of the bins are either in real-space, or based on pair-centres, neither 
of which introduces redshift distortion effects. 

We have argued, and it is clear from previous work, that any 
interpretation of projected clustering measurements must account 
for redshift space distortions. In fact, comparing correlation func- 
tions calculated using different binning schemes might actually 
prove to provide a mechanism for measuring the amplitude of the 
redshift-space distortions. This is beyond the scope of our current 
draft, and we leave this for subsequent work. 

7.2 Future Surveys 

To quantify the effect of redshift-space distortions for future sur- 
veys, we have used the expected radial selection function and pho- 
tometric redshift distribution for the Dark Energy Survey to predict 
the effect of redshift-space distortions on projected clustering mea- 
surements. This analysis is also relevant to other planned surveys 
such as PanStarrs and the LSST, which will have similar radial se- 
lection functions. We have contrasted two different types of bin- 
ning: top-hat — in which we only allow galaxies between a given 
radial bound to enter our sample — and pair-centre — in which we 
only count galaxy pairs with an average radial position that lies 
within our bounds. For typical bin widths that will be applied to 
these surveys, we find that top-hat binning in the radial direction 
leaves a strong signal from redshift-space distortions. Using a pair- 
centre binning scheme reduces the redshift-space distortion signal, 
by as much as 80% in realistic situations (see Fig. 1 14b and should 
therefore allow the measurements to be more sensitive to the cos- 
mological parameters one wishes to constrain. 

In this analysis, we have only considered the simplified situ- 
ation where the redshift-space distortions act along one axis of a 
Cartesian basis. However, the arguments we have put forward in 
favour of pair-centre binning do not rely on this assumption, and 
will remain valid even when wide-angle effects are included in any 
analysis. 
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